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Abstract 



In this paper, we derive exact closed-form expressions for the bivariate Nakagami-m cumulative 
^ distribution function (CDF) with positive integer fading severity index m in terms of a class of hy- 

, ^ pergeometric functions. Particularly, we show that the bivariate Nakagami-m CDF can be expressed 

QQ as a finite sum of elementary functions and bivariate confluent hypergeometric $3 functions. Direct 

' applications which arise from the proposed closed-form expression are the outage probability (OP) 

^ I analysis of a dual-branch selection combiner in correlated Nakagami-m fading, or the calculation of the 

level crossing rate (LCR) and average fade duration (AFD) of a sampled Nakagami-m fading envelope. 

^ Index Terms 

in 

Average fade duration, Bivariate Nakagarm-m, Correlation, Cumulative distribution function, Di- 
versity reception. Level crossing rate, Nakagami-m fading. 

o 

(N 

^ I. Introduction 

> 

•'-J There exist a number of statistical distributions which model the random fluctuations of the 

X 

^ signal amplitude when transmitted through a wireless channel [1]. Besides Rayleigh, Rician and 

other fading models, the Nakagami-m distribution [2] is an empirically-validated distribution 
which is extensively used to model different propagation conditions by means of changing the 
value of the fading severity index m. This distribution includes Rayleigh and one-sided Gaussian 
as particular cases, and approximates very closely both Rician and lognormal distributions [3]. 
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While the PDF and the CDF of a Nakagami-m variate follow very simple and tractable 
expressions, the statistical characterization of the multivariate Nakagami-m distribution is rather 
involved and remains as an open problem in the literature. Although a closed-form expression 
for the joint PDF of two correlated Nakagami-m variates was given in the original manuscript 
by Nakagami [2], a closed-form expression for the bivariate Nakagami-m CDF has not been 
obtained yet to the best of our knowledge. Existing results provide expressions for this bivariate 
CDF in terms of infinite summations [4-6], and more recently [7] in terms of a single integral 
involving special functions. The result in [7] is particularly relevant, since it allows to express 
the PDF and the CDF of a number of multivariate distributions in terms of a single integral, 
whilst previous results were only available in terms on nested infinite summations [8, 9]. 

As the Rayleigh distribution is coincident with the Nakagami-m distribution when the fading 
parameter m = 1, it may result expectable that the existing closed-form expression for the 
bivariate Rayleigh CDF [10] in terms of the first order Marcum-Q function is a particular case 
of the bivariate Nakagami-m CDF. This aspect was discussed by Simon and Alouini in [1, p. 
174], who pointed out that "One might anticipate that the bivariate Nakagami-m CDF could be 
expressed in a form analogous to ( the bivariate Rayleigh CDF) depending instead on the m}^- 

order Marcum Q-function Unfortunately, to the authors knowledge an expression analogous 

to (the bivariate Rayleigh CDF) has not been reported in the literature, and the authors have 
themselves been unable to arrive at one." . 

In this paper, we demonstrate that the bivariate Nakagami-m cumulative distribution function 
can be expressed in closed-form as a finite sum of elementary functions and $3 functions, for a 
positive integer fading severity index m. The confluent hypergeometric function of two variables 
$3, which is well studied in classical books of integrals and special functions [11, 9.261.3] and 
Laplace transforms [12], is one of the bivariate forms of Kummer's confluent hypergeometric 
function iFi and is often referred to as the Horn function [13]. Interestingly, there exists a 
connection between this kind of hypergeometric functions and the Marcum-Q function [14] of 
first-order, which allows to reduce the proposed expression for the bivariate CDF to the Rayleigh 
case when m = 1. Hence, the connection foreseen by Simon and Alouini is found at the level of 
this family of hypergeometric functions, and a further connection between $3 and generalized 
Marcum-Q functions of higher orders may be inferred. 

The bivariate confluent hypergeometric function $3 appears in the literature in a number of 
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scenarios, e.g. in the CDF of the minimum eigenvalue of a correlated non-central Wishart matrix 
[15], or in statistical models for multisensor synthetic aperture radar images [16]. Although 
$3 function is not yet included in most commercial mathematical packages, its calculation 
can be easily performed using a numerical Laplace transform inversion [1, 17]. Therefore, a 
simple Mathematica™ routine for the calculation of the confluent hypergeometric function $3 
is given. This program provides accurate results for any values of the power correlation coefficient 
< p < 1, thus circumventing some of the issues related with the convergence of the infinite 
series expressions for the CDF [4-6] as p — > 1. 

Summarizing, the main mathematical contributions of this paper are listed below 

• A closed-form expression for the bivariate Nakagami-m CDF is given, in terms of a finite 
sum of elementary functions and bivariate confluent hypergeometric functions $3. 

• The connection between the closed-form expression for the bivariate Rayleigh CDF [10] 
and the derived CDF expression with m = 1 is established. 

These relevant results find numerous applications in the field of wireless communications; 
particularly: 

(1) We provide a closed-form expression for the outage probability (OP) of a dual branch 
selection combining (SC) scheme under correlated Nakagami-m fading. 

(2) We derive closed-form expressions for the level crossing rate (LCR) and the average 
fade duration (AFD) of a sampled Nakagami-m fading envelope, following the new approach 
introduced in [18]. 

The remainder of this paper is organized as follows: in Section II, the main mathematical 
contributions of this paper are presented. Section III is dedicated to analyze some scenarios of 
interest in communications making use of the derived expressions for the Nakagami-m bivariate 
CDF. Numerical results are introduced in Section IV, whereas the main conclusions are exposed 
in Section V. 

II. Statistical Analysis 

A. Notation and preliminaries 

In the following, we use | • | to indicate the modulus of a complex number, and E{ } to 
represent the expectation operation. The correlation coefficient of two random variables X and 
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Y is defined as 

^ cov{X,Y} 

px,Y - — , ^ ^ , (i; 

where cov{} denotes co variance operation, and (j\,aY represent the variance of the random 
variables X and Y respectively. 

The CDF of a random variable X is defined as Fx{x) = Pr{X < x}, and consequently, the 
complementary CDF (CCDF) of X is defined as Fx{x) = Pr{X > x} = 1- Fx{x). 

The joint CDF of two correlated random variables X and Y is defined as 

Fx,Y{x,y)^PT{X <x,Y <y}, (2) 
whereas their joint CCDF is defined as 

FxA^^y) = MX >x,Y>y} = Fx,Y{x,y) - Fx{x) - Fyiy) + 1. (3) 

For the sake of compactness, we use the notation — Iv where Jp represents the integral 
along a contour F. 

We introduce an integral associated with the Marcum-Q function of m-th order, as well as 
some contour integrals of interest. 

Definition 1: Incomplete integral 7, (5) of Qm function. 

/•oo 

Hm{u,l,^)- / x'^'^~'^exp [-mx"^) Qm{7x,S)dx, (4) 

J u 

where Qm{a, b) is the m*'* order Marcum Q function, m e N, e [0, 00) and 5, 7 e R. 
Definition 2: Contour integral Ik 

where F is a circular contour of radius less than unity that encloses the origin, k eN,u e [0, 00) 
and (5, 7 G M. 
Definition 3: Contour integral Ji^k 

J,Au,l,5,p) 4 £^-^^p-'^exp exp dp (6) 

where F is a circular contour of radius less than unity and greater than p that encloses the origin, 

{i, k} eN,p e {0,l),u e [0, 00) and (5, 7 e R. 
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B. Results 

In this subsection, the main mathematical contributions of this paper are presented. In order 
to obtain a closed-form expression for the integral Hru', we first establish the connection with 
the integrals Xk and Ji^k- 

Proposition 1: The integral l-Lrn {u, 7, ^) can be expressed as 

^ j {m-iy.J2 («y j j^/m,k {u, S, 7) , (7) 



k=0 

„2 



where a = 7 + 2m and 

2 \ (m—k) rn—k / 2 \ (m— fe— 

\ Am, I 

i=\ 

Proof: See Appendix A. ■ 
Next, we express the integrals Xk and Ji^k in terms of the confluent hypergeometric $3 function. 

Proposition 2: The contour integral Xk can be expressed in closed-form as 



where $3 {a,b;x,y) is the confluent hypergeometric function of two variables [11, 9.261.3]. 
Proof: See Appendix B. ■ 
Proposition 3: The contour integral Ji^k can be expressed in closed-form as 

J-.. 7, P) - [j) ^s[i,k + ^■,p-,-—j, (10) 

where $3 (a, 6; x, y) is the confluent hypergeometric function of two variables. 

Proof: See Appendix C. ■ 
Note that combination of (9), (10) and (7) gives the integral Hm in closed-form in terms of 
the $3 function. 

Lemma 1: The joint CCDF of two correlated Nakagami-m random variables Ri,R2 with 
positive integer fading severity index m, power correlation coefficient p, and Qi = E{i?f}, 
Vt2 = E {i?2} power can be expressed as 



- 2m™ / ri /2mp r2 / 2m , 



Proof: See Appendix D. 
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Corollary 1: Using the results in Lemma 1 and Propositions 1-3, the joint CCDF of two 
correlated Nakagami-m random variables Ri,R2 with positive integer fading severity index m, 
power correlation coefficient p, and f^i = E{Rl}, ^2 = ^{Rl} power can be expressed in 
closed-form as 



^3 (^>^ + ^;|(i^>p(ar7fe(i^)') 



J r(k 



m—k / ox k+i—l 

mvo \ 1 



i=l 



r(A; + i) 



(12) 



Corollary 2: The joint CDF of two correlated Nakagami-m random variables Ri,R2 with 
positive integer fading severity index m, power correlation coefficient p, and Qi = E{i?f}, 
Jl2 = E power can be expressed in closed-form as 



m— 1 



exp 



Fr^,r^ (ri, r2] m, p) = 1 - 

fe=0 

- (1 -p)~''exp 



2 \ k -I 

mr^ \ 1 



1^1 y A;! 



- 1 



+ exp 



2\ k 



k\ 



m r 



mrt, 



2\ k 



^2 J T {k + 1) 



1 — p \Vti Vt-} 



^2 {i-p)'P V Vnr Vf¥ (i-p) 



m—k 



1=1 ^ ^ 



2 \ k^-%-\ 



1 



V{k^i) 



r\ r2 m 



Proof: The CDF of a Nakagami-m variate is given by 

Fx(x) 



r(m) 



(13) 



(14) 



where r(-, •) is the upper incomplete Gamma function, and r(-) is the Gamma function. In case 
m is a positive integer, (14) can be expressed as 



fe=0 ^ ^ 



(15) 
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Finally, using (15) and (12) in (3), it is obtained the closed-form expression for the bivariate 
Nakagami-m CDF in (13). ■ 

These expressions for the bivariate Nakagami-m CCDF (and thus for the CDF itself) are 
new to the best of our knowledge. These results are particularly relevant since no closed-form 
expressions were available for these functions in the literature. Expression (13) for the bivariate 
Nakagami-m CDF is given in terms of a finite sum of $3 functions, which are listed in classical 
books of integrals and special functions [11, 9.261.3], as well as in the Laplace transforms tables 
[12]. Besides, it allows for connecting the results given in this paper with the existing results 
for the bivariate Rayleigh CDF thanks to the connection between $3 and Qi functions [14], as 
follows. 

Corollary 3: The joint CDF of two correlated Nakagami-m random variables Ri,R2 with 
fading severity index m, Jli = E {R\}, = E power, and power correlation coefficient 
p, specialized at m = 1, reduces to the well-known expression for the bivariate Rayleigh CDF 
[10] as 



Fr^,r^ (ri,r2;l,p) = 1 - exp (^-^^ Qi 



r2 



exp 



Qi 



2p 



r2 



2p ri 



(16) 



where Qi{a,b) is the first-order Marcum-Q function. 

Proof: Setting m = 1 in (13), and defining u = rl/^/^^ and v = r2/\/n^ we have 



,v;l,p) = exp{-u^) - ^exp {u^ + v'^)^ 



X 



^111 f uv 

$3 1,1;:^ ,P 



2 / \ 2 

V I UV 

*3 ,P 




• (17) 



I-P \l-pj J \ I-P 

Using the relationship between the confluent hypergeometric function $3 and the first-order 
Marcum-Q function given in [14] 



Qi (a, 6) = exp - — + 



a2 fe2 



2 2 



$3 1, 1; 



a2 



2 ' 4 



(18) 



and after some straightforward manipulations, it is obtained (16). 
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III. Applications 

In this Section, we use the derived closed-form expressions for the bivariate Nakagami-m 
CDF to analyze some scenarios of interest in communications: (1) outage probability analysis 
of a dual-branch selection combining scheme in correlated Nakagami-m fading, (2) calculation 
of higher order statistics of sampled Nakagami— m fading channels. 

A. Outage probability of dual-branch SC 

Let us consider a dual-branch receiver affected by correlated Nakagami-m fading. In this 
scenario, the expression for the instantaneous SNR per branch is given by 7^ = EgRi/No, 
where i = 1,2, Eg is the transmitted symbol energy and Nq is the noise power spectral density. 
Equivalently, the average SNR per branch is given by 7j = E{R^}Es/Nq. 

In a SC scheme, the receiver chooses the branch with maximum 7^ for symbol decision 
i-e. 75C — niax{7i,72}. Hence, the OP of a dual-branch SC scheme affected by correlated 
Nakagami-m fading can be expressed in terms of the bivariate CDF as 

Pscil) = Pr{7i < 7, 72 < 7} 

= Fr,,r^ (ri = ^/7l, = ^/72; ^> p) > (19) 

where p = Pr^^r^ is the power correlation coefficient according the definition in (1), Vti = 71 
and = 72 Substituting into (13), the following closed-form expression for the OP in this 
scenario is obtained 




This closed-form expression for the OP of SC diversity in correlated non-identical Nakagami- 
m fading is equivalent to those given the literature, either in terms of infinite series [4] or in 
integral form [7]. 
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B. Higher order statistics of sampled Nakagami-m fading channels 

Second-order statistics of fading channels incorporate information related with the dynamics of 
the random processes, characterizing their evolution due to the variation along a certain dimension 
(e.g. time). Specifically, the level crossing rate (LCR) measures how often the envelope fading 
crosses a threshold value, whereas the average fade duration (AFD) informs about the amount 
of time that the envelope remains below this threshold. 

Classically, the LCR and the AFD have been calculated using Rice's approach [19], based 
upon the knowledge of the statistics of the continuous fading envelope and its time derivative. 
For the particular case of Nakagami-m fading, the higher order statistics were calculated in [20] 
in closed-form. 

Recently [18], it has been proposed an alternative formulation for the calculation of higher 
order statistics of sampled fading channels which incorporates the inherent discrete-time nature 
of fading channels due to sampling. Since the probability of missing a level crossing of the 
continuous envelope fading in a sampling interval is not negligible for a finite sampling period 
[21], the LCR for a continuous fading process (i.e., using [19]) can be seen as an upper bound 
of the associated sampled random process. Using this new approach [18], we provide exact 
closed-form expressions for the LCR and the AFD of a sampled Nakagami-m fading process. 

The LCR of a sampled random process, defined as the average rate at which the envelope R 
crosses a certain threshold u in the positive (or equivalently in the negative) direction can be 
expressed as 

where i?i = R{t), R2 = R{t + T), and T denotes the sampling period. Ri and R2 are correlated 
and identically distributed random variables, with a CDF denoted as Fji{x) = Fji-^{x) = Fr^{x), 
and joint CDF defined as Fr^^r,^{x, y) = Pr{i?i < x, i?2 < y}- 

In this scenario, the LCR of a sampled random process can be expressed in compact form as 
[18, eq. 4] 

^^^^^ ^ Fr{u) - Fr^^r^{u,u) ^ Fr{u) - Fr^^r^{u,u) ^2^^ 
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The AFD of a sampled random process, defined as the average duration of the envelope R 
remaining below a specified threshold level u, can be calculated as [18, eq. 6] 



For the particular case of Nakagami-m fading, the exact closed-form expressions for the LCR 
and the AFD can be obtained setting {ri/v^H^ = r2/v^n^ = u}^ in (13), and substituting into 
(22) and (23) respectively. These expressions for the LCR and AFD of sampled Nakagami-m 
fading are new in the literature. 



In this Section, we evaluate the expressions derived in the previous Section in some scenarios 
of interest. For the evaluation of the bivariate confluent hypergeometric function $3, we use 
the Mathematica™ program given in Appendix E. The Nakagami-m random variables used in 
Monte Carlo (MC) simulations have been generated considering that the correlation between 
the pairs of underlying Rayleigh processes composing the two Nakagami-m signals occurs only 
between in-phase components (and equivalently between quadrature components) [6]; i.e., there 
is no cross-correlation between in-phase and quadrature components. 

Fig. 1 illustrates the OP of a dual-branch SC scheme with balanced branches (i.e.,7i=72) 
affected by Nakagami-m fading, for different values of fading severity index m and power 
correlation coefficient p. The normalized SNR threshold 'jth is defined as ^th — n^'l/'Ji- It 
is appreciated how the closed-form expression and the MC simulations perfectly match. As 
expected, the OP is reduced as the fading severity index m is increased, or the correlation 
between the reception branches is reduced. 

Fig. 2 shows the OP of a dual-branch SC scheme with unbalanced branches (i.e.,7i=572) 
affected by Nakagami-m fading, for different values of fading severity index m and power 
correlation coefficient p. The normalized SNR threshold is defined as jth — "^^7/71- In this 
situation, the outage performance is lower compared to the scenario in Fig. 1. Again, theoretical 
values and MC simulations are in excellent agreement. 

'where fii = f22 = f2 and u is the threshold value normalized to v^- 
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In Fig. 3, the level crossing rate of a sampled Nakagami-m fading process is depicted, for 
different values of Doppler frequency fd and sampling period T. We have assumed that the 

underlying Rayleigh random variables have a correlation coefficient Pr{t) = Jo{27r far), which 
yields an envelope fading normalized correlation coefficient [22] given by 



where 2^1 is the Gauss hypergeometric function. It is observed that the LCR values of the 
sampled Nakagami— m envelope differ from the LCR values of the continuous Nakagami-m 
envelope [20] for low values of the threshold level. This difference is increased as the value 
of fdT product is increased, i.e., as the fading process is less oversampled. Interestingly, the 
continuous LCR is an upper bound of the actual LCR of the sampled fading process, as indicated 
in [18]. 

Finally, Fig. 4, illustrates the average fade duration of a sampled Nakagami-m fading process, 
for the same values of Doppler frequency fd and sampling period T used in Fig. 3. In this case, 
the AFD of the continuous process is a lower bound of the actual AFD. It is observed that 
the AFD of the sampled Nakagami-m process reaches an irreducible floor equal to T; this is 
in concordance with the definition of the AFD for a sampled fading process, as the minimum 
expectable duration of a fading is necessarily one sampling period interval. 



In this paper, we have derived exact closed-form expressions for the bivariate Nakagami-m 
CDF with positive integer m, in terms of the confluent hypergeometric function of two variables 
$3. These expressions are particularly convenient, since $3 function is well studied in classical 
books of integrals, special functions and Laplace transforms, and can be efficiently computed 
using a numerical Laplace inversion method. Besides, the derived CDF is shown to reduce to 
the bivariate Rayleigh CDF given in terms of the Marcum-Q function, which suggests a further 
connection between this family of hypergeometric functions and the generalized Qm functions. 

We have illustrated the applicability of the obtained expression in two scenarios of interest in 
communications: an outage probability analysis of a dual-branch SC in correlated Nakagami-m 



^ (^^)%Fi(-0.5,-0.5,m,p«)-l 




(24) 



V. Conclusion 
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fading channels, and the calculation of the higher order statistics of sampled Nakagami-m fading 
channels. Simulations corroborate the validity of the derived closed-form expressions. 
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Appendix A 
Proof of proposition 1 

We aim to find a closed-form expression for the integral Hm (7, S^u), defined as 

POO 

'Hm('W,7,5)— / x'^"^~^ exp (^—mx^) Qrn i^x, S) dx. (25) 

J u 

The Marcum-Q function can be expressed in terms of a contour integral [23] as 

Q„ (7a:, 5) = exp ^ — ) f^-;^—p'^^[-^ + ^\dp; (26) 

where F is a circular contour of radius less than unity enclosing the origin. Thus, we can express 

W„ («, 7. = exp i-m.^) {exp (- t^'j i { ^ + } *} 

(27) 

After some manipulations, and using a = 7^ + 2m, we have 
{u, 7, 5) = exp (^-^) x^-^ exp (^-^a) {/, + T^}^^ 

(28) 

Changing the integration order in (28), we obtain 

-Hm {u, 7, 5) = exp (^-^^ I jT" x^^-' exp (^-^e^ dxj ^ exp (^^^ dp- (29) 

where £ = a — ^ 
p 

The inner integral in (29) can be expressed as [24] 



f.-expf4.)^--"'2-'rf,„,^l, 
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where r(m, w) is the upper incomplete Gamma function. In case m is a positive integer, r(m, w) 
can be expressed as 

m—l u 

E'OJ 
— ; (31) 

k=0 

Substituting (30) and (31) into (29) we have 

Hm in, 7, S) = exp (m - 1)! j( cxp x 

-S)""-K-=f(-S))S(=f)"(-5)'iK» 

which can be conveniently rearranged as 

(f 2 \ / 2\ m—l ^ 2 \ 1 

-y j exp 2-^- (m - 1)! (^a^J («,5,7) 

(33) 

where 

1 1 _j. fu^l'^\ f S'^P^ 



Using partial fraction expansion, it is possible to re-express Im,k as follows, in order to separate 

the presence of the singularities at p = 1 and p = p into two different integrals. Thus 

(35) 

where 

\ (m—k) / 1 ^ m—k—i+l 



Finally, identifying the integrals (5) and (6) in (35), yields the desired expression (7). 



Appendix B 
Proof of proposition 2 

Let us express the integral Xk {u, 7, 5) in the following compact form as 

I, {u, 7, 5) = ^ F,{p) exp (^^^ dp, (37) 
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where the contour F is defined as a circular contour of radius less than unity that encloses the 
origin and 

i^i(p) = ^p-'expf^^V (38) 



p-1 V 2 P 

Interestingly, the integrand in Ik is in the form of an inverse Laplace transform. Hence, we 
aim to find a connection between the integral defined in the contour F and the general Bromwich 
integral given by 

1 pa+joD 

fig) I ^iP) (gp) dp = {F{p)} , (39) 

where F{p) is the Laplace transform of f{g). Let us consider the contour C given in Fig. 5, 
where the value of a is chosen to be at the right of all the singularities of -Fi(p). Hence, the 
contour integral along C is given by 

£ Fi (p) exp (^^^ dp=^ Fi (p) exp (^^^ + £ """^^ (^) ^"^^^ 

Using Cauchy-Goursat theorem, we can equivalently express 

£ Fi(p) exp (^^^ £ ^^^^^ (^) ^ £ ^^^^^ (^) ^"^^^ 
where Co and Ci are closed contours which enclose the singularities at p = and p = 1 
respectively. Combining (40) and (41), we can express 

1 /■"+^'°° {5'^p\ f fs'^p\ f fs'^p\ 

— J Fi (p) exp ( — 1 dp + ^ Fi (p) exp ( — 1 dp = j Fi (p) exp ( — J dp. 

c, CM ^^^^ 

The modulus of Fi(p) in is given by 

with h = M^7^/2. Using 

V-^\ ^ WJ—u ^ I' (44) 

|p-l| ||p|-l|iJ>2-R 
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and considering that 



we have 



\M = |eRc(6/p)| <g|6/K| < glVRol^ (45) 

R>Ro 



\Fi{p)\ < 2el'/^"l . (46) 

R>Ro ^ y: 

M R-l 

Thus, \Fi{p)\^^j^^je < MR~^ for some Z > on as it! — > oo, and using Jordan's lemma 
the integral equals to 0. Hence, using the residue theorem and choosing Co = F, we have 

^fc ^' '^^ = ^ fj'^ ^'^P^ (^) ~ {^'^P^ (^) } ^"^^^ 

= jr-' {F,{p)}^^^ - Res |fi(p) exp (^^) | ^ , 

where Res {F^p)}^^^ denotes the residue of F{p) at p = x. 

Expression (47) defines a relationship between the integral defined in a closed contour F 
and the inverse Laplace transform of the integrand. Using [12, 4.24.9], we have 



where $3 {a,P;x,y) is the confluent hypergeometric function of two variables [11, 9.261.3], 
which is one of the bivariate forms of the confluent hypergeometric function iFi. Using (48), 
we can express the integral 1^ as 

Appendix C 
Proof of proposition 3 

The calculation of ^ is performed using a similar procedure to that in Appendix B. Let us 
express the integral Ji^k {u,'y,d,p) in the following compact form as 

Ji,k {U, 7, ^:P)^ ^ F2,i{p) ^^^^ 
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where the contour F is defined as a circular contour of radius less than unity that encloses 
the origin and 

Again, the integrand in Ji^^ is in the form of an inverse Laplace transform. Thus, we aim 
to find a connection between the integral defined in the contour F and the general Bromwich 
integral (39). 

Let us consider the contour C given in Fig. 6, where the value of a is chosen to be at the 
right of all the singularities of F^^iij)). Hence, the contour integral along C is given by 

J F2,iip) ^^"P [~Y ) ^ 27rj J -^2,i(p) exp ( — j dp + j) F2,i{p) exp ( — 1 dp. 

c, ^^^^ 

As in (43), |-^2,i(p)|p=^eJ^ — for some Z > on C/j as — > oo, so that the integral 

equals to 0. 

Under Cauchy - Goursat theorem, we can choose a contour F as a circular contour of radius 

less than unity enclosing the origin. Since p < 1, the singularities of -p2,i(p) at p = and p = p 
are enclosed within F. Thus, we have 

1 r+^'°° , , f6^p 



F2,,(p)exp(^^^c^p = jfF2,,(p)exp(^^^dp (53) 



Therefore, the integral Ji^k can be expressed as 



Ji,k^C-'{F2,iip)}._.; (54) 

y~ 2 

The inverse Laplace transform in (54) can be calculated using [12, 4.24.9], yielding 
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Appendix D 
Proof of lemma 1 

The joint PDF of two correlated Nakagami— m variates is given in [1, 2] by 

4m (rir2)"* 



fRuR2{ri,r2\m, p) 



Im-l 



m— 1 



exp 



2 2 



.vnji^(i-p), 

where Qi — E {r?}, /n(3:^) is the n*'* order modified Bessel function of the first kind, and p is 
the power correlation coefficient between Ri and R2 according to the definition in (1). 
For the sake of compactness, we define 



n r2 



and hence this PDF can be rewritten as 



2771"^ 2m 



(i-p)' 



r(m) (vp)' 

fx,Y{x,y;m,p) = kq; (xy)"" exp {x^ + /m-i {oiyfpxy) ; 



(57) 



(58) 



The bivariate CDF is defined as 



Fx,y{u, v;m, p) ^Fr{x <u,y <v} ^ / /x,y (a;, y)ciy(ia; 

'0 Vo 



(59) 



or equivalently, the complementary CDF (CCDF) 

/•OO /'OO 

Fx,y{u, v;m, p) ^Pr{x > u,y > v} ^ / /x,y (a;, y)ciycix (60) 

J U J V 

Substituting (58) into (60), we have 

Fx,y{u, v; m, p) ^ Ka J exp j y"" exp (-^2/^) -^m-i {a^pxy) dyj 



(61) 



/"CO 



where 



T (x, = y exp -^m-i ipc^xy) dy. (62) 

Let us first find a closed-form expression for T {x,v). Performing a change of variable 

' 1 



ay"^ = 
\fady = dz 



T{x,v) = 



^ exp -— ) Im-i (tz) dz. (63) 
a I \ a \ 2 ' 
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where t — Jotpx. After some straightforward manipulations, we can express 



T(x,^) = (^)-'"-'r-^xp (^0 (^-0 exp (^-^ {tz)dz. 



(64) 

Thus, the first integral can be expressed in terms of the Marcum-Q function Qm (a, ^) i-^- 

T (x, ^) = ^ (v^^:^)"""^ exp (^^^^^ Qm (a^v^' ^v^) ■ (65) 
Substituting (65) into (61), we obtain 

Fx,y{u, v]m,p)^ k J x"^ exp i^-^^^J {y/p^V^ exp ( J Qm {xy/ap, v^/a) dx (66) 
which can be reexpressed in compact form as 

POO 

Fx,y{u, v; m,p) = K, (x/p)""^ / x^""^ exp (-mx^) Qm {^x, 6) dx, (67) 

J u 

where 

7 = .Jap, (68) 
5 = -uv^. (69) 

The integral in (67) can be identified with (4), yielding 

Fx,y{u, V] m,p) = K ( Vp)"""^ Hm {u, 7, 5) . (70) 

Finally, substituting Fr^^r^ (ri, r2; m, p) = Fx,y {u — ri/\/rt[, v — r-ij "m, p) the desired 
expression for the bivariate Nakagami— m CCDF is obtained in (11). 

Appendix E 

MATHEMATICA™ program for the computation of $3 FUNCTION 

InvLap[fO_, xO_] := Module [{f = fO, x = xO}, 
fk;[eje_] := Cos [x eje] ; 
a = 1 + 1/x; 

fg[eje_] := (2 Exp[a x])/\[Pi] f[a + I e je] ; 

psuml[l] := NIntegrate [Re [fk [eje] fg[eje]], {eje, 0, \[Pi]/(2 x) } ] ; 
psuml [i_?NumberQ] := NIntegrate [Re [fk [eje] fg[eje]], 

{eje, -(\[Pi]/(2 X)) + i \[Pi]/x , -(\[Pi]/(2 x) ) + (i + 1) \[Pi]/x}]; 
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psuml[l] + NSum[psuml [i] , {i, 1, \ [ Infinity] } , 

Method -> "AlternatingSigns " , "VerifyConvergence" -> False]] 

(* Implementation of ^3{(3,'y; zt,yt)*) 

Conf luentHypergeometricX [CapitalPhi] 3 [\ [Beta] \ [Beta]_, \ [Gamma] \ 
\[ Gamma ]_, zz_, yy_, tt_] := 
Module [{ \ [Beta] = \ [Beta] \ [Beta] , \ [Gamma] = \ [ Gamma ] \ [ Gamma ] , 
z = zz, y = yy, t = tt, f}, 

f[s_] := Gamma [\ [Gamma] ] /s'\ [Gamma] (1 - z/s) "-\ [Beta] Exp[(y)/s]; 
l/t~ (\ [Gamma] - 1) InvLap[f, t]] 
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Normalized threshold (dB) 



Fig. 1. Outage probability versus threshold level 74/1 (normalized to 71), for balanced reception 71 = 72 and different values 
of m and p. Markers correspond to Monte Carlo simulations. 
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Normalized threshold y^^ (dB) 



Fig. 2. Outage probability versus threshold level 7th (normalized to 71), for unbalanced reception 71 = 572 and different 
values of m and p. Markers correspond to Monte Carlo simulations. 
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Fig. 3. Normalized level crossing rate Nr{u)/ fd of a sampled Nakagami-m fading channel, for different values of fdT and 
m. Markers correspond to Monte Carlo simulations. 
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Normalized threshold u (dB) 



Fig. 4. Normalized average fade duration AR{u)fd of a sampled Nakagami-m fading channel, for different values of fdT and 



March 8, 2013 



DRAFT 



25 




Fig. 5. Contour integration for integral T^. 
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Fig. 6. Contour integration for integral jTi,*,. 
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